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1.  Executive  Summary 


1.1  Background 

This  ESTCP  project  evaluated  the  benefits  and  utility  of  applying  transport  optimization 
algorithms,  operable  on  desktop  computers,  versus  a  traditional  trial-&-error  approach.  The 
focus  was  on  groundwater  pump  and  treat  (P&T)  systems.  The  transport  optimization  algorithms 
link  mathematical  optimization  techniques  with  simulations  of  groundwater  flow  and 
contaminant  transport  to  determine  the  best  combination  of  well  locations  and  pumping  rates  for 
a  P&T  system.  These  mathematical  algorithms  can  potentially  contribute  to  long-term  operating 
cost  reduction  and/or  improved  performance  of  these  systems  with  respect  to  compliance 
objectives  (e.g.  achieve  cleanup  standards  in  less  time). 

Remediation  System  Evaluation  (RSE)  or  Remedial  Process  Optimization  (RPO)  provides  a 
broad  assessment  of  optimization  for  the  remedial  systems,  which  includes  system  goals  and  exit 
strategy,  below-ground  performance,  above-ground  performance,  monitoring  and  reporting,  and 
potential  for  alternate  technology.  The  pumpage  optimization  (as  demonstrated  in  this  project)  is 
a  subset  or  a  component  of  these  more  general  optimization  evaluations  for  cases  where  P&T  is 
expected  to  be  a  major  component  moving  forward.  The  pumpage  optimization  impacts 
subsurface  aspects  of  the  remedy  (e.g.,  cleanup  time,  containment)  and  in  some  cases  also 
impacts  above-ground  aspects  (treatment  plant  size/flow  rate  and  influent  concentrations),  and 
possible  monitoring  requirements  as  well. 

A  previous  project  sponsored  by  the  US  EPA  (US  EPA,  1999a,b)  demonstrated  potential 
avoidance  of  millions  of  dollars  in  O&M  costs  over  the  projected  P&T  lifetime  at  two  of  three 
sites  through  the  application  of  “hydraulic  optimization”.  Hydraulic  optimization  couples 
simpler  optimization  techniques  (linear  and  mixed-  integer  programming)  with  simulations  of 
groundwater  flow  (but  not  transport).  The  transport  optimization  techniques  that  are  the  focus  of 
this  ESTCP  project  are  potentially  more  powerful  than  the  hydraulic  optimization  techniques, 
because  they  rigorously  incorporate  predictions  of  contaminant  concentrations,  contaminant 
mass,  and/or  cleanup  duration.  However,  transport  optimization  codes  are  also  more  complex 
and  difficult  to  apply  than  hydraulic  optimization  codes. 

This  demonstration  project  was  divided  into  two  phases: 

•  Phase  1:  Pre- optimization  site  screening 

•  Phase  2:  Demonstration  of  transport  optimization  codes 

For  Phase  1,  a  spreadsheet-based  pre- optimization  screening  methodology  was  developed  and 
applied  at  eleven  existing  DoD  P&T  systems.  The  objective  of  Phase  1  was  to  provide  end-users 
with  a  framework  and  a  simple  tool  for  quickly  and  inexpensively  prioritizing  which  sites  are 
most  likely  to  benefit  from  the  application  of  transport  optimization  codes.  For  this  project, 
criteria  for  site  selection  in  addition  to  those  specified  in  screening  methodology  included  the 
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existence  of  a  flow  and  transport  model  considered  to  be  “up-to-date  and  acceptable  for  design 
purposes”  and  a  willingness  to  consider  implementing  changes  suggested  by  the  optimization 
analysis. 

For  Phase  2,  transport  optimization  was  compared  with  a  trial- &-error  approach  for  three  sites: 

•  Umatilla  Chemical  Depot,  Hermiston,  Oregon  (“Umatilla”) 

•  Tooele  Army  Depot,  Tooele,  Utah  (“Tooele”) 

•  Former  Blaine  Naval  Ammunition  Depot,  Hastings,  Nebraska  (“Blaine”) 

Both  the  Umatilla  and  Tooele  sites  have  existing  P&T  systems  in  operation,  and  the  Blaine  site  is 
a  planned  P&T  system.  The  pre- optimization  screening  methodology  developed  in  Phase  1  can 
be  used  for  both  existing  and  planned  systems,  and  the  use  of  the  optimization  algorithms  in 
Phase  2  is  also  applicable  to  both  new  and  existing  P&T  systems. 


1.2  Objectives  of  the  Demonstration 

The  primary  objective  of  this  project  is  to  demonstrate  the  cost  benefit  of  applying  transport 
optimization  codes  to  existing  P&T  systems  relative  to  the  traditional  trial- &-error  approach.  At 
each  of  the  3  sites,  the  potential  cost  savings  from  applying  transport  optimization  exceeded  the 
expected  costs  of  the  technology.  Therefore,  this  objective  has  been  met.  A  secondary  objective 
was  to  provide  each  installation  where  the  demonstration  is  performed  with  alternate  pumping 
strategies  that  would  be  feasible  and  cost-effective  to  implement.  Based  on  the  feedback  from 
the  respective  teams,  this  objective  was  partially  but  not  fully  met,  primarily  because  the 
objectives  of  the  installation  or  the  site- specific  transport  models  changed  by  the  time  the 
demonstration  was  completed.  While  the  installations  have  been  encouraged  to  implement 
optimization  suggestions  resulting  from  the  demonstration,  they  were  not  required  to  do  so  in 
order  to  participate  in  this  project.  Several  of  the  installations  have  indicated  an  interest  in 
upgrading  their  current  groundwater  models  for  potential  future  application  of  optimization 
codes  with  the  updated  models. 

It  is  noted  that  a  numerical  groundwater  model  will  never  exactly  predict  groundwater  flow  and 
contaminant  transport,  and  that  any  results  obtained  based  on  groundwater  model  predictions 
must  be  evaluated  in  that  context.  However,  those  issues  pertain  equally  to  any  design  based  on 
flow  and  transport  modeling,  whether  obtained  using  transport  optimization  algorithms  or  trial- 
&-error  techniques.  This  project  did  not  evaluate  the  impact  of  the  uncertainty  associated  with 
simulation  model  parameters  on  the  optimization  solutions.  However,  this  issue  could  be 
evaluated  in  future  projects  either  by  examining  the  impact  to  optimal  solutions  from  varying 
model  parameter  values  or  by  using  stochastic  optimization  methods  to  identify  optimal  solutions 
that  are  robust  despite  the  uncertainty.  Optimal  solutions  are  often  “at  the  edge”  of  what  is 
feasible  and  therefore  are  not  always  robust.  This  project  did  not  evaluate  the  robustness  of  the 
optimal  solutions.  One  way  to  increase  the  robustness  of  the  solutions  would  be  to  apply  a 
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“safety  factor”  to  the  optimization  problem  (i.e.,  impose  more  restrictive  constraints  than  are 
actually  required),  which  in  general  will  lead  to  more  conservative  designs. 


1.3  Regulatory  Drivers 

There  are  no  specific  regulations  that  mandate  the  use  of  optimization  methods.  However,  P&T 
systems  that  are  not  showing  sufficient  hydraulic  capture,  or  where  remedial  schedules  do  not 
seem  to  be  achievable  with  the  existing  system  design,  are  prime  candidates  for  regulatory  and 
stakeholder  interest  and  consideration  of  optimization  efforts  such  as  the  application  of 
mathematical  optimization  methods. 

Beyond  these  considerations,  there  are  no  technology- specific  regulatory  issues  that  need  to  be 
directly  addressed  beyond  those  that  constrain  the  design  and  operation  of  the  P&T  systems 
being  examined,  e.g.,  such  as  hydraulic  capture  boundaries  and  overall  revision  of  P&T  system 
objectives.  Such  regulatory  issues  were  represented  by  the  installation  and  considered  during  the 
strategic  development  of  the  mathematical  formulations  that  were  solved  using  the  transport 
optimization  algorithms.  The  ESTCP  project  team  encouraged  regulatory  participation  in  the 
process  and  for  each  demonstration  site  offered  to  help  site  personnel  communicate  with  their 
regulatory  partners  regarding  the  optimization  technology.  However,  installation  personnel  were 
ultimately  responsible  for  keeping  regulators  involved  in  the  project  to  the  extent  desirable  and 
necessary. 


1.4  Demonstration  Results 

For  all  three  sites,  there  were  two  groups  applying  optimization  algorithms  and  one  group 
applying  trial- &-error  as  a  scientific  control.  For  each  of  the  three  sites,  multiple  formulations 
were  solved  by  each  group.  In  each  and  every  case,  the  groups  applying  the  optimization 
algorithms  found  improved  solutions  relative  to  the  trial- &- error  group.  This  project  clearly 
demonstrated  that  mathematical  optimization  is  capable  of  identifying  substantially  improved 
solutions  to  real-world  problems  encountered  for  optimization  of  P&T  systems.  At  all  three 
sites,  the  potential  cost  savings  outweighed  the  expected  costs  of  applying  the  technology.  The 
solutions  found  using  transport  optimization  algorithms  were  5%  to  50%  better  than  those 
obtained  using  trial-&-error  (measured  using  optimal  objective  function  values),  with  a 
representative  improvement  of  about  20%. 

Please  note  that  optimization  results  are  not  compared  to  the  current  system.  The  reason  is  that 
the  current  system  was  not  designed  with  the  current  version  of  the  groundwater  model,  nor  was 
the  current  system  designed  to  be  optimal  for  any  of  the  formulations  solved  in  this  study. 
Therefore,  it  is  not  fair  to  compare  the  current  system  to  the  optimal  results,  and  there  are  no 
scientific  conclusions  that  can  be  gained  from  such  a  comparison.  The  focus  of  this  project,  by 
design,  is  the  comparison  of  solutions  obtained  with  transport  optimization  algorithms  versus 
trial- &-error  approach. 
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1.5  Stakeholder/End-User  Issues 


This  demonstration  was  designed  to  evaluate  the  potential  benefit  of  transport  optimization  codes 
relative  to  the  traditional  trial- &- error  approach  for  optimizing  pumping  strategies  for  a  P&T 
system.  Thus,  this  effort  primarily  examined  technical  and  cost  related  issues,  such  as  the 
conditions  necessary  for  appropriate  or  feasible  application  of  the  technology,  the  reasons  for 
differences  in  the  optimal  results  between  groups,  and  the  factors  that  influence  performance  of 
the  technology. 
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2.  Technology  Description 


2.1  Technology  Development  and  Application 

Most  P&T  systems  have  been  designed  through  the  use  of  numerical  flow  and/or  solute  transport 
simulation  models,  such  as  MODFLOW  and  MT3D.  Transport  models  (e.g.,  MT3D)  allow  for 
prediction  of  contaminant  concentrations,  contaminant  mass,  and  cleanup  times.  The  most 
common  decision  variables  determined  with  groundwater  modeling  are  well  locations,  on/off 
status  of  particular  wells,  and  the  extraction  or  injection  rates.  The  goal  is  to  identify  the  best  (or 
optimal)  combination  of  values  for  the  decision  variables,  which  traditionally  means  that  an 
“objective  function”  is  minimized  or  maximized,  and  all  “constraints”  are  satisfied.  Potential 
objective  functions  include  “minimize  total  cost”,  “minimize  cleanup  time”,  and  many  others. 
Constraints  are  limits  on  the  extraction/injection  rates,  limits  on  well  locations  or  the  number  of 
wells,  limits  on  hydraulic  heads,  limits  on  concentrations,  limits  on  total  capital  investment,  and 
many  others. 

Traditionally,  the  groundwater  simulation  model  is  run  repeatedly  to  simulate  different  pumping 
scenarios  (referred  to  as  “trial- &- error”).  Each  pumping  scenario  is  entered  manually  with 
respect  to  locations  of  wells  and  the  pumping/injection  rates.  After  the  simulation  is  completed, 
the  modeler  determines  by  inspection  if  cleanup/containment  is  achieved  and  if  all  other  design 
constraints  are  satisfied. 

This  trial-&-error  approach  relies  heavily  on  the  experience  and  insight  of  the  modeler.  Using 
this  trial- &-error  approach,  a  limited  number  of  simulations  are  performed  (typically  10  to  50) 
and  a  preferred  pumping  strategy  is  then  selected.  A  limitation  of  this  approach  is  that  only  a 
small  number  of  pumping  strategies  can  be  investigated,  and  the  objective  function  and 
constraints  are  often  not  rigorously  stated  (in  mathematical  terms).  Another  limitation  is  that  the 
nonlinear  relationship  between  pumping  rates  and  groundwater  concentrations  (i.e., 
concentrations  do  not  change  proportionally  with  pumping  rates)  makes  it  difficult  to  select 
promising  well  locations  and  pumping  rates  based  on  earlier  choices. 

Transport  optimization  codes  (the  focus  of  this  demonstration  project)  couple  transport  models 
with  nonlinear  mathematical  optimization,  to  allow  a  more  rigorous  evaluation  of  potential 
pumping  strategies  (i.e.  using  mathematical  algorithms  instead  of  manual  iteration).  Nonlinear 
optimization  algorithms  are  required  because  concentration  changes  and/or  cleanup  time  changes 
are  not  linearly  related  to  pumping  rate.  The  coupled  simulation-optimization  approach  is 
appealing  because  it  can  presumably  identify  improved  pumping  strategies  for  a  given  objective 
function  and  constraint  set  by  more  efficiently  searching  the  range  of  potential  combinations  of 
well  rates  and  locations. 

The  nonlinear  optimization  problem  that  results  from  the  transport  optimization  formulations  can 
be  conceptualized  as  a  mountain  range  with  a  series  of  peaks  and  valleys.  The  optimal  solution 
is  either  the  highest  peak  or  the  lowest  valley,  depending  on  the  nature  of  the  objective  function 


ESTCP  Final  Cost  &  Performance  Report,  11/20/03 


5 


(maximize  or  minimize).  There  are  many  algorithms  for  solving  these  nonlinear  problems. 
Traditional  approaches  use  derivatives  of  the  objective  function  and  constraints  to  go  “uphill” 
from  the  starting  point  of  the  search  until  the  peak  is  found.  These  approaches  find  only  the 
highest  peak  or  lowest  valley  nearest  the  starting  point  of  the  search.  They  are  also  difficult  to 
implement  with  complex  transport  models,  which  may  not  be  differentiable. 

A  newer  class  of  optimization  methods,  referred  to  as  “heuristic  global  optimization  methods”, 
has  emerged  in  recent  years.  These  methods  include  simulated  annealing,  genetic  algorithms, 
outer  approximation,  and  tabu  search.  These  global  methods  are  designed  to  search  the  potential 
solution  space  for  the  highest  peak  or  lowest  valley.  These  global  methods  often  require 
intensive  computational  effort,  but  have  become  more  practical  for  application  on  personal 
computers  as  computer  speeds  have  increased.  Heuristic  techniques  can  also  handle  any  form  of 
objective  function  and  constraints  and  any  type  of  simulation  model,  along  with  relatively 
straightforward  linking  of  simulation  models  with  the  optimization  algorithm.  The  transport 
optimization  codes  demonstrated  in  this  project  use  a  variety  of  heuristic  global  optimization 
methods. 

This  project  evaluates  the  utility  of  applying  transport  optimization  codes  for  optimizing 
extraction/injection  rates  and  extraction/injection  locations.  Optimization  can  potentially  result 
in  reduced  cleanup  time  and/or  reduced  life-cycle  costs.  The  objective  function  associated  with 
each  formulation  is  designed  to  be  a  metric  for  comparing  one  solution  to  another  within  the 
optimization  approach,  and  therefore  is  ideally  suited  for  measuring  performance  of  one 
pumping  strategy  versus  another  when  all  of  the  constraints  are  satisfied. 


2.2  Process  Description 

The  process  of  performing  transport  optimization  generally  consists  of  the  following  steps: 

•  Model  Development  -  Calibrate  a  groundwater  model  to  the  point  where  it  is  considered 
acceptable  for  design  purposes 

•  Develop  Optimization  Formulation(s)  -  Define  in  mathematical  terms  an  objective 
function  to  be  minimized/maximized  and  a  set  of  constraints  that  must  all  be  satisfied 

•  Solve  Optimization  Formulation(s)  -  Determine  if  there  are  any  feasible  solutions  (i.e., 
that  satisfy  all  the  constraints),  and  if  so,  determine  the  “best”  or  “optimal”  solution 


After  an  optimization  formulation  is  solved,  the  solution  often  reveals  additional  constraints  or 
changes  in  the  objective  function  that  may  better  target  the  solution  to  meet  overall  project  goals. 
Hence,  it  is  common  to  modify  one  or  more  aspects  of  the  formulation,  and  then  to  solve  the 
modified  formulation. 
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To  solve  the  optimization  formulations,  transport  optimization  codes  that  use  a  variety  of 
mathematical  algorithms  have  been  developed  (discussed  above).  These  optimization  codes  link 
to  the  transport  model  simulation  input  and  output,  and  also  require  one  or  more  input  files  that 
specify  parameters  and  options  to  be  used  by  the  optimization  code.  In  some  cases,  the  process 
of  solving  the  optimization  formulation  also  requires  the  user  to  simplify  the  optimization 
formulation  and/or  place  limits  on  the  potential  range  of  values  for  decision  variables  to  avoid 
excessive  computational  requirements.  This  might  include  limiting  the  potential  well  locations, 
only  considering  specific  pumping  or  injection  rates  (e.g.,  increments  of  50  gpm),  not  allowing 
pumping  or  injection  rates  to  change  over  time,  etc. 


2.3  Previous  Testing  of  Technology 

Since  the  1980’s,  many  researchers  have  coupled  groundwater  simulation  models  with 
mathematical  optimization  techniques  to  address  groundwater  management  issues.  Several 
universities  have  developed  transport  optimization  codes,  and  some  have  been  tested  at  actual 
field  sites.  Three  examples  of  recent  applications  of  transport  optimization  are: 

•  Utah  State  University:  Wurtsmith  Air  Force  Base 

•  University  of  Alabama:  Massachusetts  Military  Reservation  (CS-10) 

•  Utah  State  University:  Massachusetts  Military  Reservation  (CS-10) 

Each  is  described  below.  Peralta  (2001)  describes  other  recent  real-world  design  projects  using 
the  Utah  State  University  model. 

Wurtsmith  Air  Force  Base,  MI:  Optimizing  Contaminant  Mass  Removal  Using  Artificial  Neural 
Network  -  Utah  State  University. 

In  this  case  (, Aly ,  A.H.  and  R.C.  Peralta,  1997),  transport  optimization  was  used  to  develop  an 
optimal  strategy  for  remediating  TCE  and  DCE  groundwater  plumes.  Management  goals  and 
restrictions  were  identified  and  prioritized  as  follows: 

•  Capture  the  TCE  and  DCE  dissolved  phase  groundwater  plumes 

•  Reduce  TCE  and  DCE  concentrations  to  less  than  94  ppb  and  less  than  230  ppb, 
respectively,  within  6  years 

•  Total  extraction  of  groundwater  cannot  exceed  400  gpm 

•  No  treated  water  may  be  injected  into  the  groundwater 

•  Treatment  facility  effluent  cannot  exceed  5  ppb  of  TCE 
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An  artificial  neural  network  was  used  to  simulate  contaminant  concentrations  in  the  optimization 
model.  The  model  considered  a  total  of  24  potential  extraction  well  locations.  Six  alternative 
optimal  pumping  strategies  were  ultimately  evaluated  for  the  final  design.  After  discussions  with 
stakeholders,  a  final  strategy  was  chosen  based  on  its  minimization  of  total  pumping  rates, 
minimization  of  total  time  to  meet  objectives,  and  overall  benefit  to  the  stakeholders. 

Chemical  Spill- 10  (CS-10)  site  located  at  the  Massachusetts  Military  Reservation  -  University  of 
Alabama  and  Utah  State  University 

Two  of  the  three  recent  study  applications  of  transport  optimization  were  applied  for  the  CS-10 
plume  at  Massachusetts  Military  Reservation.  A  P&T  system  is  operating  to  remediate  and 
contain  a  TCE  plume  approximately  17,000  feet  long,  6,000  feet  wide,  and  up  to  140  feet  thick. 
Between  Fall  1999  and  Spring  2000,  transport  optimization  codes  were  utilized  to  maximize 
TCE  mass  removal  over  a  30- year  time  horizon,  subject  to  the  following  constraints:  (1)  the  TCE 
concentration  must  be  lower  than  or  equal  to  5  ppb  beyond  the  base  boundary,  (2)  all  extracted 
water  must  be  reinjected  into  infiltration  trenches,  (3)  individual  wells  are  subject  to  pumping 
capacities,  and  (4)  the  total  pumping  rate  should  be  restricted  for  cost  considerations.  The 
decision  variables  were  the  extraction  rates  and  well  locations  for  four  perimeter  wells  that  were 
being  considered,  and  the  extraction  rates  for  five  in-plume  wells  that  were  already  constructed. 

Results  for  the  two  optimization  studies  are  summarized  below: 

•  University  of  Alabama.  In  this  case  (AFC  EE.  1999;  Zheng  and  Wang,  2002),  the  optimal 
strategy,  as  determined  by  the  simulation- optimization  analyses,  suggests  using  only  one 
perimeter  well  (rather  than  four  wells)  and  a  maximum  pumping  rate  of  2700  gallons  per 
minute  (gpm).  The  results  of  the  analysis  demonstrate  that  it  is  possible  to  remove  more 
TCE  mass  (approximately  3.5%)  under  the  same  amount  of  pumping  assumed  in  the  trial- 
&-error  design,  and  that  it  can  also  lead  to  substantial  cost  savings  by  reducing  the 
number  of  wells  needed  and  adapting  dynamic  pumping.  Preliminary  cost  estimates 
indicated  that  this  strategy  would  yield  life  cycle  cost  savings  of  $2.4  million.  Some 
elements  of  the  design  were  implemented. 

•  Utah  State  University.  In  this  case  ( Peralta  et  al,  1999a,  b),  the  simulation- optimization 
modeling  enhanced  mass  removal  rates  and  aided  in  well  placement,  with  an  additional 
constraint  of  preventing  the  plume  from  contaminating  clean  aquifer  between  the  western 
and  central  lobes.  Specifically,  the  modeling  identified  a  configuration  that  would  extract 
approximately  6%  more  mass  over  30  years,  while  reducing  the  extraction  rate  by  50 
gpm  and  could  cost  $0.54  million  less  in  construction  cost  alone.  With  slight  tweaking, 
this  design  was  constructed  and  is  functioning  as  expected. 
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2.4  Advantages  and  Limitations  of  the  Technology 

A  properly  defined  optimization  problem  can  be  solved  through  manual  trial- &-error  adjustment 
or  using  a  formal  optimization  technique.  While  the  trial-&-error  method  is  simple  and  widely 
used,  it  is  usually  limited  in  practice  to  a  small  number  of  simulations  (typically  10-50)  because 
it  is  labor  intensive.  The  transport  optimization  codes  more  efficiently  search  the  potential 
solution  space,  such  that  thousands  of  simulations  are  typically  performed  automatically,  and 
each  successive  round  of  new  simulations  is  designed  to  be  “more  promising”  than  the  previous 
round. 

Key  advantages  of  transport  optimization  codes  include  the  following: 

•  many  more  combinations  of  extraction  and  injection  well  rates  can  be  evaluated  using 
search  algorithms  that  are  far  more  efficient  than  trial- &- error  or  random  search 

•  the  process  of  mathematically  specifying  an  objective  function  and  a  set  of  constraints  is 
required  for  transport  optimization,  and  this  process  (frequently  overlooked  during  trial- 
&- error  modeling)  forces  competing  goals  and  strategies  to  be  considered  and  compared 

•  because  it  is  more  automated  than  trial- &- error,  transport  optimization  is  less  prone  to 
bias  in  selecting  well  rates  and  well  locations,  and  is  therefore  more  likely  to  discover 
unexpected  solutions 

Limitations  of  transport  optimization  codes  include: 

•  the  site  must  develop  a  transport  model  that  is  considered  a  reasonable  predictor  for 
design  purposes  (also  true  for  trial- &- error) 

•  the  complexity  of  applying  the  nonlinear  transport  algorithms  may  require  specialized 
expertise  for  many  real-world  groundwater  modeling  problems 

•  the  codes  are  very  computer  intensive  for  most  real-world  groundwater  modeling 
problems,  potentially  requiring  simplification  of  the  simulation  model  and/or  dedicated 
use  of  one  or  more  computers;  these  simplifications  require  specialized  expertise 

A  limitation  that  pertains  to  both  trial- &-error  and  the  use  of  transport  optimization  algorithms  is 
that  the  optimal  results  are  based  on  model  predictions,  which  are  subject  to  uncertainty.  A 
number  of  approaches  exist  for  considering  uncertainty  in  the  optimization  process,  but  these 
were  not  evaluated  in  this  project. 
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3.  Demonstration  Design 


3.1  Performance  Objectives 

The  objective  of  the  overall  project  is  to  demonstrate  the  cost  benefit  of  applying  transport 
optimization  codes  by  addressing  the  following  questions: 

•  Do  the  results  obtained  from  these  optimization  software  packages  (e.g.  recommended 
optimal  P&T  scenarios)  differ  substantially  from  the  optimal  solutions  determined  by 
traditional  “trial- &- error”  optimization  methods? 

•  Do  the  results  obtained  from  these  optimization  software  packages  warrant  the  additional 
effort  and  costs  when  compared  to  traditional  “trial- &- error”  optimization  methods? 

Three  optimization  formulations  were  developed  for  each  of  the  three  sites,  based  on  interaction 
with  the  installation.  Two  modeling  groups  used  their  own  independently  developed  transport 
optimization  software,  and  the  other  group  used  a  traditional  “trial- &- error”  optimization 
method.  The  results  from  two  separate  transport  optimization  software  programs  can  be 
compared  to  each  other,  and  to  the  trial-&-error  group,  to  assess  performance  objectives. 

Table  3- 1  summarizes  performance  criteria  discussed  in  the  Demonstration  Plan,  including  the 
expected  performance  to  be  achieved. 


Table  3-1.  Performance  Objectives  for  Overall  Project 


Type  of 
Performance 
Objective 

Primary  Performance 
Criteria* 

Expected  Performance* 
(Metric) 

Performance 

Objective 

Met? 

Qualitative 

1 .  Reduce  annual  operating  costs 

Annual  operating  costs  are 
reduced 

** 

2.  Faster  remediation 

Increased  contaminant  removal 
efficiency 

** 

3.  Reduce  cost  of  system  life 
cycle 

Reduced  annual  cost  and/or 
reduced  cleanup  time 

Yes 

Quantitative 

1 .  Reduce  annual  operating  costs 

>  20%**** 

** 

2.  Reduce  system  life  cycle  costs 

>  20%**** 

*Based  on  a  comparison  of  results  obtained  with  optimization  algorithms  versus  trial-&-error 


**As  discussed  in  section  4.3,  the  cases  based  on  cost  objectives  were  formulated  in  terms  of  life-cycle  costs, 
which  incorporate  tradeoffs  between  annual  costs  and  cleanup  time.  Therefore,  only  life -cycle  cost  reduction 
was  directly  evaluated.  In  some  cases  life-cycle  costs  reductions  were  due  to  lower  annual  costs  (e.g.,  Blaine), 
and  in  some  cases  life-cycle  cost  reductions  were  due  to  reductions  in  cleanup  time  (e.g.,  Umatilla) 
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***Metric  achieved  for  multiple  formulations,  but  not  for  all  formulations 

****  The  criterion  of  >20%  was  originally  in  the  Demonstration  Plan  Part  II.  It  was  based  on  general  experience 
in  applying  optimization  to  groundwater  modeling  problem  provided  by  discussions  with  GeoTrans,  Dr. 

Peralta,  and  Dr.  Zheng. 


The  last  column  of  Table  3- 1  indicates  whether  or  not  the  performance  objectives  were  met 
according  to  project  results. 


3.2  Selection  of  Test  Sites 

In  this  project,  a  screening  method  was  developed  for  site  selection  (see  Appendix  G  of  the 
Technical  Summary  Report).  The  screening  analysis  is  a  two-stage  procedure.  The  first  stage  is 
intended  to  quickly  remove  sites  from  consideration  if  they  are  not  likely  to  benefit  from  either 
hydraulic  or  transport  optimization,  based  on  the  following  simple  questions: 

•  Are  O&M  costs  >  $100K/year? 

•  Is  the  system  flowrate  >  50gpm? 

•  Is  the  estimated  cleanup  time  >  5  years? 

If  the  answers  to  all  three  questions  are  “Yes”,  a  potential  benefit  from  hydraulic  and/or  transport 
optimization  is  suggested,  and  the  second  stage  (i.e.,  quantitative  potential  cost  saving 
evaluation)  is  recommended  to  classify  the  sites  into  Tiers  regarding  potential  benefits  that  might 
be  realized  by  performing  a  hydraulic  optimization  or  a  transport  optimization  analysis. 

Three  sites  were  ultimately  selected  for  the  transport  optimization  demonstration: 

•  Umatilla  Chemical  Depot,  Hermiston,  Oregon  (“Umatilla”) 

•  Tooele  Army  Depot,  Tooele,  Utah  (“Tooele”) 

•  Former  Blaine  Naval  Ammunition  Depot,  Hastings,  Nebraska  (“Blaine”) 

Umatilla  and  Tooele  have  existing  P&T  systems  in  operation,  and  Blaine  is  in  the  design  stage 
for  a  planned  P&T  system.  Table  3-2  summarizes  the  results  of  the  screening  analysis  for  these 
three  sites. 


Table  3-2.  Selection  Criteria  of  Demonstrations  Sites 


Selection  Criteria 

Umatilla 

Tooele 

Blaine 

O&M  Costs  >  $100K/yr? 

$430K/yr 

$lM/yr 

$2M/yr 

System  Flowrate  >  50  gpm? 

1,300  gpm 

>  5,000  gpm 

4,000  -  5,000  gpm 

Estimated  Cleanup  Time  >  5  Yrs? 

10-30  years 

>  5  years 

50  -  80  years 
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Max.  Potential  Cost  Savings  (Life- 
Cycle)  -  Hydraulic  Optimization  * 

$415,485 

$3,379,423 

$11,488,043 

Max.  Potential  Cost  Savings  (Life- 
Cycle)  -  Transport  Optimization  * 

$362,986  - 
$1,298,000 

$3,329,423  - 
$5,683,687 

$11,435,543  - 
$14,359,313 

*  These  are  “pre-optimization  ”  estimates,  not  actual  optimization  results. 

*  For  calculations  of  potential  cost  savings,  it  was  assumed  that  the  system  durations  of  hydraulic  optimization 
for  Umatilla,  Tooele,  and  Blaine  were  2  0  years,  20  years,  and  50  years,  respectively.  For  the  transport 


optimization,  calculations  were  based  on  a  10-30%  reduction  in  cleanup  duration. 


For  the  purpose  of  this  demonstration,  two  other  criteria  were  used  for  site  selection,  which  were 
met  by  all  three  sites: 

•  Up-to-date  flow  and  transport  models  exist  that  are  considered  reasonable  to  apply  for 
design  purposes  at  the  site 

•  Site  managers  expressed  a  willingness  to  consider  implementing  the  recommendations 
that  might  arise  from  the  optimization  results 


3.3  Test  Site  History/Characteristics 

Brief  descriptions  of  the  test  facilities  are  provided  below.  More  detail  is  provided  in  the 
Technical  Summary  Report. 

3.3.1.  Umatilla  Chemical  Depot  (“Umatilla”) 

Umatilla  is  located  in  northeastern  Oregon,  three  miles  south  of  the  Columbia  River  and  six 
miles  west  of  Hermiston,  Oregon.  From  the  1950s  until  1965,  the  depot  operated  an  onsite 
explosives  washout  plant.  The  plant  processed  munitions  to  remove  and  recover  explosives 
using  a  pressurized  hot  water  system.  The  wash  water  from  the  plant  was  disposed  in  two 
unlined  lagoons,  located  northwest  of  the  plant,  where  wash  water  infiltrated  into  the  soil. 
Explosives  contained  in  the  wash  water  migrated  into  the  soil  and  groundwater  at  the  site.  Two 
of  the  most  common  contaminants  in  groundwater  are  RDX  and  TNT.  These  constituents  are 
used  as  indicator  parameters  because  they  are  found  at  high  concentrations  relative  to  other 
parameters.  Remediation  of  the  groundwater  is  scheduled  to  continue  until  the  concentratio  n  of 
explosives  in  the  aquifer  meets  cleanup  levels.  The  cleanup  level  for  RDX  is  2.1pg/l  and  for 
TNT  is  2.8  (lg/1. 

The  hydrogeology  for  Umatilla  consists  of  an  alluvial  aquifer  overlying  silt  and  weathered  basalt. 
The  RDX  plume  is  much  bigger  in  area  than  the  TNT  plume.  This  is  because  TNT  is  more 
strongly  absorbed  to  the  aquifer  materials,  and  therefore  its  movement  is  retarded  relative  to 
groundwater  velocity  to  a  much  greater  extent  than  RDX.  A  groundwater  P&T  system  was 
implemented  in  January  1997.  Design  of  the  groundwater  treatment  system  was  based  in  part  on 


ESTCP  Final  Cost  &  Performance  Report,  11/20/03 


12 


the  results  of  groundwater  modeling  studies.  The  current  P&T  system  has  three  active  extraction 
wells  and  three  active  infiltration  basins.  The  infiltration  basins  are  located  around  the  perimeter 
of  the  pre-pumping  RDX  plume,  and  were  intended  to  augment  hydraulic  control.  An  additional 
infiltration  basin  (the  location  of  a  former  industrial  lagoon)  was  used  as  an  infiltration  basin  in 
the  early  stages  of  the  remedy  to  promote  in-situ  flushing  of  the  unsaturated  zone.  However,  use 
of  that  location  for  infiltration  was  terminated  due  to  concerns  that  it  could  cause  unwanted 
spreading  of  the  TNT  plume.  The  contaminated  groundwater  is  extracted  from  the  wells  and  then 
sent  to  GAC  units,  which  remove  the  contaminants.  The  treated  water  is  discharged  to  the 
infiltration  basins.  The  annual  O&M  cost  for  the  current  system  is  approximately  $430K/year. 

3.3.2.  Tooele  Army  Depot  (“Tooele”) 

Tooele,  located  several  miles  south  of  the  Great  Salt  Lake  in  Utah  was  established  in  1942  to 
provide  storage,  maintenance  and  demilitarization  of  troop  support  equipment,  especially 
wheeled  vehicles  and  conventional  weapons.  From  1942- 1966,  large  quantities  of  hazardous 
materials  were  used  and  generated  from  these  operations  in  the  industrial  area,  and  discharged 
into  ditches  and  ponds.  Ultimately  groundwater  was  impacted,  and  the  primary  contaminant  of 
concern  is  TCE,  which  was  used  as  a  solvent  in  the  repair  operations  of  military  equipment. 

The  aquifer  of  concern  generally  consists  of  alluvial  deposits.  However,  there  is  an  uplifted 
bedrock  block  at  the  site  where  groundwater  is  forced  to  flow  from  the  alluvial  deposits  into 
fractured  and  weathered  rock  (bedrock),  and  then  back  into  alluvial  deposits.  The  uplifted 
bedrock  block  and  adjoining  low  hydraulic  conductivity  alluvium  are  the  hydraulically 
controlling  features  of  the  study  area  due  to  the  steep  gradients  they  cause.  Flow  through  the 
bedrock  block  consists  of  a  steep  gradient  when  entering  the  bedrock,  a  flatter  gradient  through 
the  bedrock  core  and  a  steep  gradient  when  exiting  the  bedrock.  The  impacts  have  been  divided 
into  a  “Main  Plume”  and  a  “Northeast  Plume”.  This  optimization  study  pertained  to  the  Main 
Plume.  Concentrations  are  significantly  lower  in  the  deeper  portions  of  the  aquifer  than  in 
shallow  portions  of  the  aquifer.  Also,  the  extents  of  the  shallow  and  deep  plumes  do  not  directly 
align,  indicating  a  complex  pattern  of  contaminant  sources  and  groundwater  flow.  Continuing 
sources  of  dissolved  contamination  are  believed  to  exist. 

A  P&T  system  has  been  operating  since  1993.  The  system  consists  of  16  extraction  wells  (15 
operating  and  one  not  operating)  and  13  injection  wells.  An  air-stripping  plant,  located  in  the 
center  of  the  plume,  is  capable  of  treating  8000  gpm  of  water  (currently  treats  about  5000  gpm). 
Based  on  the  well  locations  and  previous  plume  delineations,  the  original  design  was  for  cleanup, 
but  subsequently  it  was  determined  that  the  source  area  extended  further  to  the  south.  As  a  result, 
the  current  system  essentially  functions  as  a  containment  system  (there  are  no  extraction  wells  in 
the  area  of  greatest  contaminant  concentration).  Historically,  the  target  containment  zone  has 
been  defined  by  the  5  ppb  TCE  contour.  Given  the  current  well  locations  and  continuing 
sources,  anticipated  cleanup  time  is  “a  very  long  time”.  However,  a  revised  (i.e.,  smaller)  target 
containment  zone  is  now  being  considered,  based  on  risks  to  potential  receptors.  A  revised  target 
containment  zone  might  correspond  to  the  20  ppb  or  50  ppb  TCE  contour.  Annual  O&M  cost  is 
approximately  $1M. 


ESTCP  Final  Cost  &  Performance  Report,  11/20/03 


13 


3.3.3.  Former  Blaine  Naval  Ammunition  Depot  (“Blaine”) 

Blaine  is  located  immediately  east  of  Hastings,  Nebraska,  and  was  built  in  the  early  1940s  as  an 
active  “load,  assemble,  and  pack”  ammunition  facility  during  World  War  II  and  the  Korean 
Conflict.  Waste  materials  were  generated  through  discharge  of  wastewater  to  surface 
impoundments  and  natural  drainage  areas  of  the  facility,  and  disposal  of  solid  waste  and 
explosives.  Beginning  in  the  mid-  1960s,  large  tracts  of  the  former  depot  were  either  sold  to 
various  individuals,  businesses,  and  municipalities  or  transferred  to  other  governmental  agencies. 
With  sale  and  transfer  of  the  land  to  the  U.S.  Department  of  Agriculture  (USDA)  and  area 
farmers,  over  100  irrigation  wells  have  been  installed  on  the  former  depot.  Groundwater 
contamination  at  Blaine  was  discovered  in  the  mid-1980s.  The  remedial  investigation  and  the 
annual  groundwater  sampling  results  identified  seven  source  areas  for  VOCs  with  plumes 
commingling  at  six  of  the  source  areas  and  one  primary  source  for  explosives.  Extensive 
remediation  of  source  areas  by  soil  vapor  extraction  (SVE)  or  soil  excavation  is  being 
implemented  or  has  been  completed. 

Groundwater  is  encountered  in  the  study  area  approximately  100  feet  below  ground  surface.  The 
three  saturated  hydrogeologic  units  of  primary  interest  of  this  study  are  an  unconfined  aquifer,  a 
confining  layer,  and  a  semi- confined  aquifer.  During  irrigation  season,  which  lasts  about  two 
and  half  months,  heavy  pumping  from  extensive  irrigation  wells  dramatically  alters  the 
groundwater  flow  direction.  VOC  plumes  encompass  nearly  six  and  one- half  square  miles. 
Additionally,  groundwater  contamination  from  explosives  extends  over  an  area  of  approximately 
three  square  miles  and  is  commingled  with  the  VOC  plume(s)  in  several  areas. 

There  is  no  existing  groundwater  extraction  remediation  system  at  Blaine.  This  site  is  in  the 
conceptual  design  stages,  based  on  a  draft  Feasibility  Study  performed  in  August  2000.  The 
optimization  project  is  restricted  to  simulation  of  two  parameters.  Site  managers  selected  TCE 
and  TNT  as  the  parameters  most  important  to  remedial  design.  However,  site  managers  also 
indicated  a  preference  to  not  ignore  the  other  parameters.  Therefore,  an  approach  was  developed 
to  indirectly  incorporate  the  distribution  of  other  constituents  of  concern. 


3.4  Setup  and  Operation 

The  traditional  trial- &- error  method  was  used  by  GeoTrans  to  serve  as  a  scientific  control  for  the 
transport  optimization  groups.  Two  transport  optimization  modeling  groups,  Dr.  Chunmiao 
Zheng  of  University  of  Alabama  (UA)  and  Dr.  Richard  Peralta  of  Utah  State  University  (USU), 
used  their  own  independently  developed  simulation-optimization  software  for  this  study.  These 
investigators  were  chosen  based  upon  existence  of  optimization  packages  and  prior  field 
implementation  of  their  optimization  packages  in  a  similar  fashion  to  what  was  intended  for  this 
project,  although  the  specific  codes/algorithms  they  would  apply  for  this  project  were  their 
choice. 
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Once  sites  were  selected  (Phase  1),  the  following  activities  were  associated  with  performing  the 
transport  optimization  for  each  site  (Phase  2): 

•  An  initial  draft  of  potential  optimization  formulations  was  developed  by  the  ESTCP 
project  team,  based  on  a  site  visit  (performed  in  Phase  1)  and  subsequent  phone 
conversations  and/or  e-  mails  with  the  installation 

•  Feedback  on  the  initial  draft  optimization  formulations  was  provided  by  the  installation, 
including  details  on  cost  coefficients  and/or  constraint  values 

•  The  simulation  models  were  modified  (as  necessary)  to  require  no  more  than  2  hours  of 
computational  time  and  to  include  no  more  than  2  constituents 

•  The  formulations  were  finalized  and  distributed  to  each  modeling  group  by  GeoTrans, 
including  a  feasible  solution  if  one  had  been  determined  during  the  formulation  process;  a 
FORTRAN  post-processor  for  determining  the  objective  function  value  and  status  of  the 
constraints  for  any  specific  combination  of  well  rates  simulated  with  the  transport  model 
was  also  provided  by  GeoTrans 

•  Optimization  of  each  of  the  three  formulations  for  the  site  was  performed  over  a  period  of 
approximately  four  months,  during  which  time  the  three  modeling  groups  were  not 
allowed  to  discuss  their  progress  with  each  other  or  with  the  installation  (bi-weekly 
progress  reports  were  submitted  by  each  group  to  the  US  Navy) 

•  After  the  optimization  period  for  a  specific  site,  each  modeling  group  submitted  a  report 
describing  the  results  for  each  formulation 

•  The  ESTCP  project  team  met  to  present  and  interpret  results,  with  a  subsequent 
presentation  of  results  to  the  installation  by  a  subset  of  the  ESTCP  project  team 

A  summary  of  these  activities  for  each  of  the  three  demonstration  sites,  including  the  schedule,  is 
presented  in  Table  3-3. 


Table  3-3.  Summary  of  Activities 


Activity 

Umatilla 

Tooele 

Blaine 

Site  Visit  (Phase  1) 

8/23/00 

5/31/01 

10/17/01 

Develop  Formulations 

10/16/00-3/21/01 

6/4/01  -  10/31/01 

1/15/02  -  5/15/02 

Optimization  Period 

3/22/01  -  7/16/01 

11/1/01  -  2/28/02 

5/17/02  -  9/17/02 

Project  Team  Meeting  to  Present 
Results 

10/18/01 

3/20/02 

9/19/02 
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Present  Results  to  Installation 

11/15/01 

5/16/02 

9/30/02 

Follow-up  with  Installation 

3/19/02 

12/02 

12/02 

A  third-party  expert  (Dr.  Barbara  Minsker)  was  added  to  the  project  team  during  the  beginning 
of  Phase  2  to  evaluate  and  interpret  the  results  and  their  meaning  relative  to  overall  project 
objectives. 

Brief  descriptions  of  the  setup  of  the  optimization  problems/formulations  are  provided  below. 
More  detail  is  provided  in  the  Technical  Summary  Report. 


Umatilla 


A  simple  description  of  the  formulations  is  as  follows: 


Formulation  1 :  Minimize  the  life-cycle  cost  (until  cleanup  of  both  RDX  and  TNT), 
subject  to:  1)  the  current  capacity  of  the  treatment  plant  is  held 
constant,  2)  the  cleanup  of  both  RDX  and  TNT  is  within  20  years. 

Formulation  2 :  same  as  Formulation  1,  but  allows  the  capacity  of  the  treatment  plant  to 

increase  to  a  maximum  of  1950  gpm. 


Formulation  3 :  Minimize  the  total  mass  remaining  (RDX  plus  TNT)  in  layer  1  within 
20  years. 


Tooele 


A  simple  description  of  the  formulations  is  as  follows: 


Formulation  1 :  Minimize  total  cost  over  21  years,  subject  to  meeting  specific 

concentration  limits  at  the  POE  (Point  of  Exposure)  boundary,  which  is 
located  along  a  portion  of  the  downgradient  property  boundary. 

Formulation  2 :  Same  as  formulation  1,  but  also  meet  aggressive  concentration  limits  at 
POC  (Point  of  Compliance)  boundaries,  which  are  located  in  the 
interior  portion  of  the  site. 


Formulation  3 :  Same  as  Formulation  2,  but  include  a  declining  source  term  rather  than 
a  continuous  source  term  for  unremediated  sources,  and  add  additional 
cleanup  constraints  within  the  plume. 
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Blaine 


A  simple  description  of  the  formulations  is  as  follows: 


Formulation  1 :  Minimize  life-cycle  cost  (until  cleanup)  subject  to:  1)  the  plumes 

cannot  spread  above  cleanup  levels  beyond  specified  areas;  2)  cleanup 
of  both  TCE  and  TNT  must  be  within  30  years  in  model  layers  3-6. 

Formulation  2 :  same  as  Formulation  1  but  assumes  diversion  of  2400  gpm  of  extracted 
water  to  a  utility  (i.e.,  that  water  does  not  require  treatment  and 
subsequent  discharge). 


Formulation  3 :  Minimize  the  maximum  total  remediation  pumping  rate  in  any 

management  period  over  a  30-  year  simulation,  such  that  the  plumes  do 
not  spread  above  cleanup  levels  beyond  specified  areas. 


A  brief  summary  of  the  formulations  for  each  site  is  provided  in  Table  3-4.  The  detailed  formulations  are 
provided  in  the  Appendices  of  the  Technical  Summary  Report. 


Table  34.  Formulation  Summary  (Key  Aspects)  for  the  Three  Demonstration  Sites 


Site  Name 

Objective  Function 

Major  Constraints 

Umatilla 

Form.  1 

Minimize  life-cycle  cost 
until  cleanup 

1.  Current  treatment  capacity 

2.  Cleanup  of  RDX  and  TNT  <  20  yrs 

Form.  2 

Minimize  life-cycle  cost 
until  cleanup 

1.  Increased  treatment  capacity 

2.  Cleanup  of  RDX  and  TNT  <  20  yrs 

Form.  3 

Minimize  total  mass 
remaining  in  layer  1  after 
20  yrs 

1 .  Cleanup  of  RDX  and  TNT 

2.  Limit  on  #  new  wells  and  recharge  basins 

Tooele 

Form.  1 

Minimize  total  cost  over 

21  years 

1.  POE  concentration  limit  for  TCE  at  site 
boundary  after  3  yrs 

Form.  2 

Minimize  total  cost  over 

21  years 

1.  POE  concentration  limit  for  TCE  at  site 
boundary  after  3  yrs 

2.  POC  concentrations  limits  for  TCE  at  specific 
locations/times  within  site  boundary 

Form.  3 

Minimize  total  cost  over 

21  years 

1.  POE/POC  concentration  limits 

2.  Declining  source  term 

3.  Cleanup  (TCE  <  50ppb)  at  most  locations  <  9  yrs 

Blaine 

Form.  1 

Minimize  life-cycle  cost 
until  cleanup 

1.  Plume  containment 

2.  Cleanup  of  TCE  and  TNT  <  30  yrs 

Form.  2 

Min  life-cycle  cost  until 
cleanup  w /  2400  gpm 
extracted  water  diverted 

1.  Plume  containment 

2.  Cleanup  of  TCE  and  TNT  <  30  yrs 

Form.  3 

Minimize  maximum  total 
pumping 

1.  Plume  containment 

2.  Limit  on  number  of  new  wells 

*Note:  see  Appendices  D-F  in  Technical  Summary  Report  for  detailed  formulations  for  each  site 
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4.  Performance  Assessment 


4.1  Performance  Data 

To  assess  performance  of  the  transport  optimization  algorithms  relative  to  the  conventional  trial- 
&-error  approach,  the  optimization  formulations  were  solved  by  the  UA  team  and  the  USU  team 
using  optimization  algorithms,  and  by  the  GeoTrans  team  using  trial-&-error  as  a  scientific 
control.  The  results  were  then  compared.  The  detailed  results  comparison  and  individual 
modeling  group  reports  are  included  in  the  Technical  Summary  Report.  The  comparative  results 
are  summarized  below. 

Umatilla  Performance  Data 


Table  4-1  shows  the  results  for  all  three  formulations.  The  objective  for  Formulations  1  &  2  is  to 
minimize  the  life-cycle  cost  until  cleanup  of  both  RDX  and  TNT.  The  objective  for  Formulation 
3  is  to  minimize  mass  remaining  (RDX  plus  TNT)  in  model  layer  1  after  20  years.  The  detailed 
formulation  document  is  included  in  the  Technical  Summary  Report. 


Table  4-1.  Umatilla  Results  Summary 


Transport  Optimization 
Algorithms 

Trial-& -Error 

Percentage 
Improvement  * 

UA 

USU 

GeoTrans 

Formulation  1: 

Objective  Function  Value 

$1.66M 

$1.66M 

$2.23M 

23% 

Formulation  2: 

Objective  Function  Value 

$1.66M 

$1.66M 

$2.02M 

15% 

Formulation  3: 

Objective  Function  Value 

0.19  kg 

0.20  kg 

0.38  kg 

50% 

*  Percentage  improvement  is  for  transport  optimization  algorithms  compared  to  trial-&-error 


Tooele  Performance  Data 


Table  4-2  shows  the  results  for  all  three  formulations.  The  objective  for  all  three  formulations  is 
to  minimize  the  total  cost  over  21  years.  The  USU  team  did  not  submit  a  design  for  Formulation 
2  as  posed  because  they  added  a  constraint  to  prevent  potential  for  mass  migration  around  one 
concentration  constraint  boundary  defined  in  the  formulation.  Therefore,  the  USU  results  for 
Formulation  2  cannot  be  directly  compared  to  the  results  from  the  other  groups.  All  three  teams 
reported  that  Formulation  3  was  infeasible  as  stated,  due  to  the  restriction  on  the  number  of  new 
wells  that  could  be  installed.  The  detailed  formulation  document  is  included  in  the  Technical 
Summary  Report. 
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Table  4-2.  Tooele  Results  Summary 


Transport  Optimization 
Algorithms 

Trial-&-Error 

Percentage 
Improvement  * 

UA 

USU 

GeoTrans 

Formulation  1: 

Objective  Function  Value 

$12.67M 

$14.14M 

$14.63M 

3%  -  13% 
(11%  -42%)** 

Formulation  2: 

Objective  Function  Value 

$14.45M 

*** 

$16.32M 

11% 

(30%)** 

Formulation  3 : 

Objective  Function  Value 

NA 

*P ercentage  improvement  is  for  transport  optimization  algorithms  compared  to  trial-&-error 


**P ercentage  in  parentheses  is  calculated  after  removing  ~  $10M  fixed  cost  that  could  not  change  with  pumping 
strategy  due  to  the  fixed  system  duration,  The  $10M  fixed  cost  includes  ~$7M  of  fixed  O&M  cost  and  ~$3M  of 
the  sampling  cost  that  cannot  be  reduced  based  on  feedback  from  the  optimization  modelers 
***  The  USU  team  did  not  submit  a  design  for  Formulation  2  as  posed  because  they  added  a  constraint  to  prevent 
potential  for  mass  migration  around  the  west  side  of  POC-MP1,  so  their  results  cannot  be  compared  directly  to 
the  other  groups 

****N0  solutions  could  be  found  that  satisfied  all  the  constraints 

Blaine  Performance  Data 


Table  4-3  shows  the  results  for  all  three  formulations.  The  objective  for  Formulations  1  &  2  is  to 
minimize  the  life- cycle  cost  until  the  cleanup  of  both  TCE  and  TNT  within  30  years.  The 
objective  for  Formulation  3  is  to  minimize  the  maximum  total  remediation  pumping  rate  in  any 
management  period  over  a  30- year  simulation.  The  detailed  formulation  document  is  included  in 
the  Technical  Summary  Report. 


Table  4-3.  Blaine  Results  Summary 


Transport  Optimization 
Algorithms 

Trial-& -Error 

Percentage 
Improvement  * 

UA 

USU 

GeoTrans 

Formulation  1: 

Objective  Function  Value 

$45.28M 

S40.82M 

S50.34M 

10%  -  20% 

Formulation  2: 

Objective  Function  Value 

S24.04M 

$18.88M 

S28.39M 

15%  -  33% 

Formulation  3 : 

Objective  Function  Value 

2737  gpm 

2139  gpm 

2879  gpm 

5%  -  26% 

*  Percentage  improvement  is  for  transport  optimization  algorithms  compared  to  trial-&-error 
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4.2  Performance  Criteria 


Primary  and  secondary  performance  criteria  for  the  optimization  analyses  are  presented  in  Table 
4-4.  These  criteria  were  applied  to  the  results  from  each  of  the  optimization  analyses. 


Table  44.  Performance  Criteria 


Performance 

Criteria 

Description 

Primary  or 
Secondary 

Reduce  Annual 
Operating  Costs 

Does  demonstration  indicate  potential  for  reducing  annual 
operating  costs  (based  on  modeling)? 

Primary 

Faster 

Remediation 

Does  demonstration  indicate  potential  for  a  shorter  duration  of 
pump  and  treat  operations  to  accomplish  a  comparable  level  of 
cleanup  (based  on  modeling)? 

Primary 

Reduce  Life  Cycle 
Cost  of  System 

Does  demonstration  indicate  potential  for  reduced  life  cycle 
based  on  capital  costs,  modified  annual  costs,  and  modified 
operating  (based  on  modeling)? 

Primary 

Factors  Affecting 
Technology 
Performance 

Extent  to  which  site -specific  factors  affect  technology 
performance  (or  prohibit  application  of  the  technology),  such  as 
reliability  of  models,  confidence  in  plume  delineation,  confidence 
in  source  area  delineation,  etc. 

Secondary 

Ease  of  Use 

What  is  the  required  skill  level  and  training  required  to  apply  the 
technology  at  other  sites,  and  can  others  be  expected  to  apply 
technology  as  effectively  (and  for  similar  cost)  as  the  project 
team  for  this  demonstration  project? 

Secondary 

Primary  performance  criteria  were  assessed  based  on  values  of  the  objective  functions  for 
competing  solutions.  The  secondary  performance  criteria  are  qualitative  rather  than  quantitative. 
Table  4-5  shows  the  primary  and  secondary  performance  criteria  along  with  expected 
performance  metrics  and  performance  confirmation  methods. 


Table  4-5.  Expected  Performance  and  Performance  Confirmation  Methods 


Performance  Criteria 

Expected 

Performance 

Metric 

Performance  Confirmation  Method 

Primary  Criteria  (Quantitative) 

Reduce  Annual  Operating  Costs 

>20% 

Objective  function  and/or  constraint  set 
evaluated  using  available  groundwater  model 

Faster  Remediation 

>20% 

Objective  function  and/or  constraint  set 
evaluated  using  available  groundwater  model 
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Reduce  Life  Cycle  Cost  Of 
System 

>20% 

Objective  function  and/or  constraint  set 
evaluated  using  available  groundwater  model 

Secondary  Criteria  (Qualitative) 

Factors  Affecting  Performance 

No  metrics 
assumed 

Feedback  obtained  from  three  demonstration 
sites 

Ease  of  Use 

Useful  to  transport 
simulation 
modelers 

Experience  of  the  project  team  from 
application  of  codes  at  three  demonstration 
sites 

Post-modeling  adjustments  are  beyond  the  scope  of  this  study  and  installations  are  not  required 
to  implement  modifications  based  on  the  demonstration  project  results.  Therefore,  the 
performance  evaluation  in  this  section  relies  solely  upon  the  most  currently  available  models  at 
the  time  of  the  implementation  effort  and  optimization  results  for  this  effort,  and  not  data  from 
future- planned  or  since  completed  adjustments  to  the  groundwater  model  and/or  to  the  pump  and 
treat  system  design. 

Please  note  that  optimization  results  are  not  compared  to  the  current  system.  The  reason  is  that 
the  current  system  was  not  designed  with  the  current  version  of  the  groundwater  model,  nor  was 
the  current  system  designed  to  be  optimal  for  any  of  the  formulations  solved  in  this  study. 
Therefore,  it  is  not  fair  to  compare  the  current  system  to  the  optimal  results,  and  there  are  no 
scientific  conclusions  that  can  be  gained  from  such  a  comparison.  The  focus  of  this  project,  by 
design,  is  the  comparison  of  solutions  obtained  with  transport  optimization  algorithms  versus 
trial-  &- error. 


4.3  Data  Assessment 

In  each  and  every  case,  the  groups  applying  the  optimization  algorithms  found  improved 
solutions  relative  to  the  trial-&-error  group.  Because  multiple  sites  were  evaluated,  and  multiple 
formulations  for  each  site  were  evaluated,  there  is  a  high  degree  of  confidence  in  the  conclusion 
that  the  application  of  optimization  algorithms  provides  improved  solutions  for  problems  posed 
in  the  manner  demonstrated  in  this  project  (i.e.,  mathematical  formulations  with  an  objective 
function  to  be  minimized/maximized  and  a  series  of  constraints).  Each  of  the  five  performance 
criteria  listed  in  Section  4.2  are  evaluated  below. 

Reduce  Annual  Operating  Costs  (Quantitative) 

During  the  formulation  process  for  each  site,  minimizing  annual  O&M  cost  was  discussed  but 
never  selected.  At  each  site  minimizing  life-cycle  costs  rather  than  O&M  costs  was  always 
preferred,  because  life- cycle  costs  can  potentially  be  reduced  even  if  annual  O&M  costs  increase 
(by  reducing  the  cleanup  time).  Evaluating  the  life-cycle  cost  considers  the  tradeoff  between  up¬ 
front  costs,  annual  O&M  costs,  and  the  cleanup  time.  Therefore,  the  performance  criterion 
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relating  to  life-cycle  costs  (discussed  below)  was  determined  to  be  more  applicable.  The  results 
indicated  that  life- cycle  cost  could  be  minimized  at  Umatilla  by  minimizing  the  cleanup  duration 
rather  than  minimizing  annual  costs,  at  Tooele  by  minimizing  the  number  of  new  extraction 
and/or  injection  wells  installed  (i.e.,  capital  costs)  rather  than  minimizing  annual  costs/pumping 
rates,  and  at  Blaine  by  minimizing  the  annual  O&M  costs  rather  then  shortening  the  cleanup 
time. 

Faster  Remediation  (Quantitative) 

During  the  formulation  process,  minimizing  cleanup  time  was  discussed  for  the  Umatilla  and 
Blaine  sites  but  never  selected.  At  each  site  minimizing  life-cycle  costs  rather  than  cleanup  time 
was  always  preferred,  because  life- cycle  costs  can  potentially  be  reduced  even  if  cleanup  time 
increases  (by  reducing  annual  or  up-front  costs).  Evaluating  the  life-cycle  cost  considers  the 
tradeoff  between  up-front  costs,  annual  O&M  costs,  and  the  cleanup  time.  Therefore,  the 
performance  criterion  relating  to  life-cycle  costs  (discussed  below)  was  determined  to  be  more 
applicable. 

However,  minimizing  life-cycle  costs  within  Umatilla  Formulations  1  and  2  also  effectively 
minimized  the  cleanup  time.  For  Formulation  1  the  transport  optimization  algorithms  identified 
solutions  with  cleanup  time  in  4  years,  versus  6  years  for  trial- &-error  (a  33%  improvement).  For 
Formulation  2,  all  three  groups  obtained  the  solutions  cleanup  in  4  years,  but  the  transport 
optimization  algorithms  achieved  lower  life-cycle  costs  than  the  trial-&-error  group.  Also,  in 
Formulation  3  for  Umatilla  (Table  4-1),  the  mass  remaining  in  solutions  from  transport 
optimization  groups  is  approximately  50%  less  than  that  of  the  trial- &-error  group.  This  also 
represents  faster  remediation.  These  results  indicate  the  potential  for  transport  optimization  to 
provide  greater  than  20%  faster  remediation  (the  metric  in  Table  4-5).  For  Blaine,  all  three 
groups  found  that,  although  reduction  in  cleanup  time  was  possible,  the  least  cost  solutions  came 
from  minimizing  pumping  in  each  management  period  rather  than  further  shortening  the  cleanup 
duration.  This  supports  the  conclusion  that  reducing  life-  cycle  cost  is  a  more  general 
performance  criterion  than  minimizing  annual  costs  or  cleanup  time. 

Reduce  Fife  Cycle  Cost  of  System  (Quantitative) 

In  most  of  the  formulations,  the  objective  was  to  minimize  life-cycle  costs.  As  listed  in  Tables  4- 
1  to  4-3,  the  transport  optimization  algorithms  frequently  determined  solutions  with  more  than 
20%  life-cycle  cost  reduction  relative  to  trial-&-error  (the  metric  in  Table  4-5),  and  20%  appears 
to  be  a  representative  value  achieved  for  life-cycle  cost  reduction  relative  to  trial-&-error  for 
these  cases  (the  representative  savings  are  somewhat  higher  than  20%  if  fixed  costs  for  the 
Tooele  site  are  excluded  from  the  formulation,  as  indicated  on  Table  4-2).  The  differences  in  the 
results  between  the  two  optimization  groups  is  likely  due  to  one  or  more  of  the  following  factors 
(detailed  in  the  Technical  Summary  Report): 

•  Different  approaches  taken  to  overcoming  the  computational  barriers  of  solving  these 
complex  problems; 
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•  Simplifications  that  individual  modelers  made  in  the  formulations  (primarily  additional 
constraints)  to  overcome  perceived  problems  in  the  solutions  they  obtained. 

It  is  considered  less  likely  by  the  project  team  that  the  differences  are  due  to  convergence  of  the 
heuristic  optimization  algorithms  to  sub -optimal  solutions 

Factors  Affecting  Performance  (Qualitative) 

After  obtaining  the  optimization  results  from  this  demonstration  project,  the  site  managers  of  all 
three  sites  chose  to  improve  the  underlying  flow  and  transport  model  and/or  further  delineate  the 
plumes  or  source  areas  prior  to  implementing  optimization  solutions.  This  highlights  that,  like 
trial- &- error  optimization,  these  algorithms  are  only  as  good  as  the  underlying  model.  In  all 
cases,  there  were  simplifications  made  in  the  formulation  process.  For  instance,  the  cost  of  a 
new  well  was  approximated  without  location- specific  details  such  as  exact  well  depth  or  piping 
costs  included.  This  highlights  that  some  simplifications  are  required  in  formulating  the 
optimization  problems,  and  therefore  many  different  alternative  formulations  can  be  developed 
(e.g.,  with  different  simplifications,  different  cost  coefficients,  etc.).  Also,  the  demonstration 
project  required  the  modeling  groups  to  perform  the  optimization  over  a  fixed  period  of  time, 
with  no  contact  with  the  installations.  In  many  cases,  questions  arose  from  the  initial  solutions 
developed  with  the  optimization  algorithms,  and  the  modelers  expressed  that  they  would  have 
preferred  to  iterate  to  an  improved  formulation  based  on  contact  with  the  installation.  Therefore, 
the  project  demonstrates  the  value  of  such  iterations. 

Ease  of  Use  (Qualitative) 

This  project  demonstrated  that  applying  the  transport  optimization  algorithms  to  these  complex, 
real-world  sites  was  more  than  just  “hitting  the  go  button”.  Primarily,  it  required  expertise  to 
limit  the  potential  solution  space  to  be  searched.  The  transport  optimization  teams  employed 
sequential  solution  approaches  to  reduce  computational  effort,  in  which  some  parts  of  the 
problem  were  fixed  while  others  were  optimized.  In  some  cases,  problems  were  solved  one 
management  period  at  a  time,  and/or  determining  well  locations  by  first  assuming  steady- state 
pumping  rates  and  then  optimizing  well  rates  over  time  for  those  pre- determined  well  locations. 
These  approaches  require  substantial  expertise  and  professional  insight. 


4.4  Technology  Comparison 

The  results  clearly  indicate  that  mathematical  optimization  methods  are  able  to  identify  solutions 
that  are  better  than  those  obtained  using  traditional  trial-  &-error  approaches.  The  solutions  found 
were  3%  to  50%  better  than  those  obtained  using  trial-&-error  (measured  using  optimal  objective 
function  values),  with  a  representative  improvement  of  about  20%.  Given  the  computational 
limits,  it  is  not  practical  for  the  optimization  algorithms  to  search  the  entire  solution  space  so 
there  is  always  a  chance  that  global  optimal  solution  will  be  missed,  but  for  the  problems  solved 
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in  this  project  better  solutions  were  consistently  found  by  the  transport  optimization  algorithms 
than  with  trial- &- error. 

This  project  did  not  specifically  address  optimization  results  versus  the  current  system,  because 
the  current  systems  were  not  designed  with  the  current  groundwater  model  and/or  were  not 
designed  based  on  the  same  optimization  formulations  considered  in  this  demonstration  project 
(i.e.,  different  objective  functions  and  constraints).  For  example,  the  pump  and  treat  systems 
might  have  been  installed  to  achieve  hydraulic  capture  of  plume  migration  whereas  the 
demonstration  project  may  have  focused  on  cleanup;  therefore,  the  optimization  results  would 
not  be  directly  comparable  to  the  existing  design  for  plume  containment.  Therefore  emphasis  on 
such  comparisons  is  not  appropriate.  However,  if  such  comparisons  were  to  be  made  by 
calculating  the  objective  function  value  (using  the  optimization  formulations  for  the 
demonstration  project)  based  on  the  current  systems  (e.g.,  existing  well  rates  and  well  locations), 
improved  objective  values  for  the  optimization  results  would  be  evident.  For  instance,  for 
Umatilla  Formulation  #1,  the  objective  function  value  for  the  current  system  (e.g.,  existing  well 
rates  and  well  locations)  would  be  $3.8M,  versus  $1.66M  from  the  optimization  results.  This  is 
largely  due  to  the  improved  cleanup  time  (17  years  for  the  current  system  versus  4  years  for  the 
optimization  results).  However,  the  trial- &- error  group  achieved  a  solution  of  6  years  for 
cleanup  time,  with  life-cycle  cost  of  $2.23M,  indicating  that  much  of  the  improvement  associated 
with  the  optimization  results  relative  to  the  current  system  is  not  attributed  solely  to  the 
optimization  algorithms.  In  this  example,  the  benefit  of  the  optimization  algorithm  (versus  trial- 
&-error)  was  to  lower  the  simulated  cleanup  time  from  6  years  to  4  years,  thereby  lowering  life- 
cycle  costs  from  $2.23M  to  $1.66M. 

The  challenges  in  applying  optimization  algorithms  increased  with  the  complexity  of  the  site 
hydrogeological  features  and  contamination.  The  greatest  challenge  the  optimization  modeling 
teams  faced  was  the  computational  requirements  of  the  optimization  algorithms  and  underlying 
simulation  models.  If  a  single  optimization  run  were  set  up  to  solve  the  entire  problem  as 
formulated,  with  all  possible  pumping  rates  and  well  locations  in  all  potential  management 
periods,  the  number  of  decision  variables  would  be  much  larger  and  the  computational  times 
associated  with  the  optimization  algorithms  would  be  prohibitive  on  today’s  computers.  Instead, 
the  teams  limited  potential  well  locations  and/or  employed  sequential  solution  approaches  to 
reduce  computational  effort,  in  which  some  parts  of  the  problem  were  fixed  while  others  were 
optimized.  These  approaches  require  substantial  expertise  and  professional  insight. 

A  limitation  of  the  trial-&-error  approach  is  that  the  objectives  and  constraints  are  often  not 
rigorously  stated.  Another  limitation  is  that  there  are  an  infinite  number  of  combinations  for  well 
locations  and  well  rates  that  are  possible,  but  the  trial-&-error  method  is  practically  limited  to 
only  a  small  number  of  numerical  simulations  (typically  10  -  50).  The  transport  optimization 
codes  more  efficiently  evaluate  the  potential  solution  space,  such  that  many  more  combinations 
(i.e.,  thousands)  of  extraction  and  injection  well  rates  and  locations  can  be  evaluated.  Also, 
transport  optimization  is  less  prone  to  bias  in  selecting  well  rates  and  well  locations  because  it  is 
more  automated  than  trial-  &- error,  and  therefore  is  more  likely  to  discover  unexpected  solutions. 
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This  project  did  not  evaluate  the  impact  of  uncertainty  in  model  parameter  values  on  the  results 
of  the  optimization  solutions.  However,  this  issue  could  be  evaluated  in  future  projects  either  by 
examining  the  impact  to  optimal  solutions  from  varying  model  parameter  values  or  by  using 
stochastic  optimization  methods  to  identify  optimal  solutions  that  are  robust  despite  the 
uncertainty. 
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5.  Cost  Assessment 


5.1  Cost  Reporting 

This  project  demonstrated  the  cost  benefit  of  applying  computer-based  transport  optimization 
codes,  and  the  cost  reporting  for  such  an  effort  differs  from  other  ESTCP  cleanup  category 
technology  demonstration  projects.  Some  typical  cost  tracking  categories  do  not  apply,  while 
there  are  other  costs  unique  to  a  modeling  optimization  demonstration  effort.  The  majority  of  the 
costs  related  to  this  demonstration  were  labor  costs  of  the  modelers.  There  are  no  capital  costs 
associated  with  this  demonstration  since  the  optimization  codes  and  existing  models  run  on 
standard  PCs.  The  project  cost  is  summarized  in  Table  5-1. 


Table  5-1.  Project  Cost  Summary  For  Contractors 


Costs  Associated  With  Three  Demonstration  Sites 

Umatilla 

Tooele 

Blaine 

Sub-Total 

Pre -Optimization  Tasks 

1)  Site  Visit  and/or  Transfer 
Information 

GeoTrans 

$3,500 

$3,400 

$3,800 

$10,700 

UA 

$13,796 

$14,526 

$14,526 

$48,848 

USU 

$9,141 

$9,664 

$9,664 

$28,469 

2)  Model  Modifications  * 

$7,730 

$7,000 

$7,000 

$21,730 

3)  Develop  3  Optimization  Formulations 

$16,000 

$14,000 

$12,000 

$42,000 

Optimization  Modeling 

4)  Solve  Optimization 

Formulations 

GeoTrans 

$32,000 

$30,000 

$31,500 

$93,500 

UA 

$75,996 

$60,345 

$60,345 

$196,686 

USU 

$56,152 

$52,467 

$53,427 

$162,046 

Reporting 

5)  Prepare  Report  and/or  Present 
Results 

GeoTrans 

$10,000 

$10,500 

$9,000 

$29,500 

UA 

$24,602 

$24,602 

$24,602 

$73,806 

USU 

$9,426 

$9,756 

$9,756 

$28,938 

*Model  modifications  were  performed  by  the  UA  group  to  for  the  demonstration  purpose  of  this  project 
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These  contractor  costs  form  the  basis  for  determining  what  future  applications  of  the  technology 
might  cost  (Section  5.2).  To  execute  this  project,  additional  costs  were  incurred  for  project 
management  and  support,  including  the  following  activities: 


•  Preparation  of  Tech  Demo  Plans  for  Phase  1  and  Phase  2 

•  Site  selection  activities  and  the  development  of  the  site  screening  methodology  referred 
to  as  “Phase  1”  of  the  project 

•  Site  visits  to  candidate  sites  that  were  ultimately  not  selected  (George,  Comhusker,  and 
Shaw) 

•  Meetings  with  installations  to  present  the  optimization  results 

•  Project  management  and  oversight,  including  contracting,  project  meetings  and 
conference  calls 

•  Preparation  of  ESTCP  Technical  Report  and  Cost- Performance  Reporting 

•  Briefings  to  ESTCP 

•  Technology  transfer  activities 

These  project  management  costs  for  executing  the  entire  project  are  summarized  below: 


Table  5  -2.  Cost  for  Project  Management,  Entire  Project 


Agencies 

Costs 

Navy  -  NFESC 

$170,000* 

EPA  -  TIO 

$27,000 

USACE  -  HTRW-CX 

$39,500  (by  8/31/03) 

USACE- WES 

$45,000 

GeoTrans,  Inc 

$129,300** 

Dr.  Minsker 

$40,000 

USACE  Site  Specific  Project  Teams 

Estimate  $17,500 

*  Cost  includes  NFESC  labor,  travel,  other  direct  costs,  and  contract  admin  charges  ( estimated  through 
completion  of  project).  As  of  8/31/03  the  amount  was  $142K,  estimating  another  $28K  to  completion 
**  As  of  8/31/03  the  total  GeoTrans  cost  was  $257,000,  estimating  another  $48K  to  project  completion  for  a  total 
of  $305,000.  Removing  the  items  from  Table  5-1  ($175,700),  the  GeoTrans  costs  related  to  project  management 
are  estimated  at  $129,300. 


5.2  Cost  Analysis 

The  major  driver  for  the  transport  optimization  modeling  is  the  model  execution  time  of  one  flow 
and  transport  simulation  for  the  underlying  model,  which  includes  model  computational  time  for 
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simulating  a  chemical  constituent  and  the  number  of  chemical  constituents  that  must  be 
simulated  to  adequately  address  the  plume  management  issues  at  the  site.  Table  5-3  lists  the 
chemical  constituents  simulated,  computational  time  for  each  flow  and  transport  simulation  for 
each  site  (not  including  computational  time  for  optimization),  and  approximate  number  of 
optimization  simulations  performed  for  each  formulation.  Generally  “simulations”  refers  to  the 
number  of  iterations  of  the  groundwater  model.  However,  due  to  the  use  of  substituted  functions 
in  place  of  the  numerical  model  in  some  formulations,  it  is  impossible  to  exactly  calculate  the 
number  of  completed  groundwater  model  simulations  performed  for  the  optimization  codes. 


Table  5-3.  Computational  Time  And  Approximate  Number  of  Simulations 


Constituents 

Simulated 

Computational 
Time  For  Each 
Flow  &  Transport 
Simulation 

Number  of  Simulations  Performed  For  Each 
Formulation 

Transport  Optimization 
Algorithms 

UA  and  USU  Teams 

Trial-& -Error 

GeoTrans 

Umatilla 

2 

(RDX  &  TNT) 

~  10  minutes 

~  1000  -  8000 
simulations 

-25-40  simulations 

Tooele 

1 

(TCE) 

~  10  minutes 

Up  to  8000  simulations 

-  60  -  80  simulations 

Blaine 

2 

(TCE  &  TNT) 

~  2  hours 

~  Hundreds/thousands 
simulations* 

-  60  simulations 

UA  group  used  less  accurate  but  much  faster  solvers  in  MT3DMS  in  the  early  stage  so  that  many  model 
simulations  were  much  shorter  than  2  hours 


As  illustrated  in  Table  5-3,  several  thousand  flow  and  transport  simulations  are  usually 
performed  using  the  transport  optimization  algorithms  to  achieve  optimal  or  near- optimal 
solutions.  Therefore,  the  computational  time  for  performing  one  flow  and  transport  simulation  is 
critical  during  transport  optimization  modeling.  Based  on  this  project,  it  is  recommended  that  the 
underlying  flow  and  transport  model  be  modified  to  the  extent  possible  prior  to  optimization  to 
shorten  the  computational  time  without  losing  model  accuracy. 

Based  on  the  competitive  bids  evaluated  in  this  project  for  selecting  the  transport  optimization 
groups,  plus  the  costs  associated  with  GeoTrans’  participation  in  the  project,  the  expected  costs 
(and  expected  time  duration)  of  applying  this  technology  at  a  future  site  is  approximated  in  Table 
5-4.  The  estimated  range  in  costs  results  from  differing  site  and  model  complexities.  The  costs  to 
conduct  transport  optimization  in  Table  5-4  are  primarily  governed  by  the  simulation  time  and 
not  any  specific  aspect  of  site’s  hydrogeologic  complexity. 
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Table  54.  Approximate  Cost  To  Apply  Transport  Optimization  Algorithms  at  A  Site 


Costs  Associated  With  Basic  Items  * 

Low  Cost 

Typical  Cost 

High  Cost 

Expected 

Duration 

Al)  Site  Visit  and/or  Transfer 
Information 

$2,500 

$5,000 

$10,000 

1-2  months 

A2)  Develop  3  Optimization 
Formulations 

$5,000 

$10,000 

$15,000 

1-2  months 

A3)  Solve  Optimization 
Formulations 

$25,000 

$40,000 

$60,000 

2-4  months 

A4)  Prepare  Report  and/or 
Present  Results 

$5,000 

$15,000 

$25,000 

1  month 

A5)  Project  Management 
and/or  Administration 

$2,500 

$5,000 

$10,000 

NA 

Total 

$40,000 

$75,000 

$120,000 

5-9  months 

Costs  Associated  With  Optional  Items 

Low  Cost 

Typical  Cost 

High  Cost 

Expected 

Duration 

Bl)  Update  and  Improve 

Models 

0 

$20,000 

$50,000 

Add  1-3  months 

B  2)  Up  To  3  Additional 
Formulations 

$15,000 

$25,000 

$40,000 

Add  2-3  months 

B3)  Additional  Contaminant 

$10,000 

$20,000 

$30,000 

Add  1-2  months 

B4)  Transport  Simulation  of  3 
hrs  each  (i.e.,  1  hr  longer) 

$10,000 

$20,000 

$30,000 

Add  1-2  months 

*  assumes  1  or  2  constituents,  and  simulation  time  of  2  hours  or  less 


Note  that  the  actual  cost  spent  by  the  two  transport  optimization  groups  solving  three 
formulations  in  the  demonstration  project  (Item  4  in  Table  5-1)  was  generally  $50K  to  $60K 
(and  in  one  case  more  than  $75K),  which  is  at  the  high  estimated  cost  for  solving  three 
optimization  formulations  (Item  A3  in  Table  5-4).  This  is  likely  attributable  to  the  fact  that  some 
cost  was  incurred  by  the  optimization  groups  to  improve  their  transport  optimization  codes  while 
solving  the  optimization  problems  during  the  demonstration  project,  and  also  because  of  the 
competitive  nature  of  the  demonstration  project,  which  probably  led  to  more  effort  trying  to 
achieve  a  global  optimum  than  would  normally  be  expended. 

With  respect  to  potential  cost  savings  to  the  government,  assume  700  P&T  systems  at  Superfund 
sites  ( USEPA ,  2003),  with  an  estimated  annual  O&M  cost  of  $570K/yr  ( USEPA ,  2003),  yielding 
a  total  O&M  cost  of  $400M/yr  (700systems  x  $570K/system).  If  the  typical  site  has  a 
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remediation  timeframe  of  20  years,  that  represents  lifecycle  costs  of  approximately  $8  billion, 
non- discounted  ($400M/yr  *  20  yrs).  If  it  is  assumed  that  10%  of  systems  have  transport 
optimization  potential,  and  if  it  is  further  assumed  that  transport  optimization  might  save  20%  of 
the  life- cycle  costs  at  those  sites,  then  the  potential  lifecycle  cost  savings  to  the  government  can 
be  estimated  to  be  $160M,  non- discounted  (i.e.,  $8  billion  *  0.1  *  0.2).  The  costs  of 
implementing  this  type  of  effort,  as  outlined  in  Table  5-4,  are  not  included  within  these  cost 
savings  projections. 


5.3  Cost  Comparison 

The  appropriate  way  to  compare  the  costs  of  applying  transport  optimization  algorithms  with 
conventional  trial- &-error  in  this  demonstration  project  is  to  compare  the  cost  associated  with 
Item  A3  in  Table  5-1,  which  is  “solving  optimization  formulations”  or  searching  for  the  optimum 
solution  based  on  mathematically  stated  formulations  (objective  functions  and  constraints).  The 
cost  comparison  associated  with  solving  optimization  formulations  and  the  percentage 
improvement  in  objective  functions  for  demonstration  sites  are  listed  in  Table  5-5.  Two  teams 
using  transport  optimization  algorithms  spent  from  approximate  $50,000  to  $76,000  to  solve 
three  optimization  formulations,  which  is  compared  to  the  actual  costs  of  approximately  $32,000 
for  the  trial- &-error  effort  for  each  of  the  three  sites  examined  on  this  project  (approximately  the 
same  amount  was  budgeted  for  trial- &-error  at  each  site,  so  for  more  complex  models  such  as 
Blaine  fewer  trial- &- error  simulations  were  performed). 


Table  5-5.  Cost  Comparison  vs.  Percentage  Improvement  For  Demonstration  Sites 


Costs  Associated  With  Solving  Optimization 
Formulations 

Percentage  Improvement  in 
Objective  Functions 

Transport  Optimization  Algorithms  vs. 
Trial  & -Error 

Triaf&-Error 

Transport  Optimization 
Algorithms 

Umatilla 

$32,000 

$  75,996  (UA) 

$  56,152  (USU) 

23%,  15%,  50%  (UA) 

23%,  15%,  47%  (USU) 

Tooele 

$30,000 

$  60,345  (UA) 

$  52,467  (USU) 

13%,  11%,  -  /(42%,  30%,  -)*  (UA) 
3%,  -,-/(ll%,  -,  -)*  (USU) 

Blaine 

$31,500 

$  60,345  (UA) 

$  53,427  (USU) 

10%,  15%,  5%  (UA) 

19%,  33%,  26%  (USU) 

**Percentage  in  parentheses  is  calculated  after  removing  ~  $  10M  fixed  cost  that  could  not  change  with  pumping 
strategy  due  to  the  fixed  system  duration,  The  $10M  fixed  cost  includes  ~$7M  of  fixed  O&M  cost  and  ~$3M  of 
the  sampling  cost  that  cannot  be  reduced  based  on  feedback  from  the  optimization  modelers 


The  university  team  costs  in  Table  5-5  likely  included  some  code  development  that  would  not 
typically  be  a  part  of  an  optimization  application.  Hence,  we  expect  that  a  comparison  of 
expected  costs,  shown  in  Table  5-4,  to  benefits  is  more  appropriate  for  estimating  future  costs 
and  benefits  of  the  technology. 
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An  issue  is  the  extent  to  which  application  of  transport  optimization  algorithms  cost  more  than 
the  application  of  trial- &- error,  and  to  compare  that  with  the  anticipated  benefits  that  might  be 
afforded  by  the  application  of  transport  optimization  algorithms  versus  use  of  trial- &- error. 

As  shown  in  Table  5-4,  the  estimated  cost  of  applying  transport  optimization  algorithm  (Item 
A3)  for  problems  like  those  formulated  for  this  project  is  approximately  $25,000  to  $60,000  (i.e., 
up  to  2  constituents,  simulations  up  to  2  hrs  long,  up  to  3  formulations).  The  cost  for  the  trial- &- 
error  group  for  Item  A3  for  this  project  was  approximately  $30,000  per  site,  although  that  group 
reported  for  each  site  that  it  would  have  performed  fewer  simulations  if  not  done  within  the 
context  of  this  demonstration  project.  Thus,  it  is  assumed  that  for  comparable  projects  (i.e.,  up  to 
2  constituents,  simulations  up  to  2  hrs  long,  up  to  3  formulations)  trial- &-error  may  cost 
approximately  $20,000  to  $25,000.  Therefore,  the  premium  for  applying  the  transport 
optimization  may  be  as  little  as  zero,  or  as  much  as  $40,000. 

The  improvements  in  objective  function  values  achieved  in  this  demonstration  project  with 
transport  optimization  algorithms  (versus  trial-&-error  methodology)  range  from  3%  to  50% 
with  the  fixed  cost  associated  with  the  Tooele  site  included,  which  can  not  change  with  the 
pumping  strategy  due  to  the  fixed  system  duration,  and  range  (for  all  three  sites)  from  5%  to  50% 
with  that  fixed  cost  removed  from  the  calculation,  with  a  typical  value  of  20%  to  30%. 

Assuming  the  objective  function  is  in  terms  of  cost,  the  potential  life-cycle  cost  savings 
associated  with  the  application  of  transport  optimization  algorithms  will  almost  certainly  exceed 
the  premium  of  up  to  $40,000  for  applying  the  technology  at  most  sites  that  satisfy  the  simple 
site- screening  criteria  (more  than  $  100,000/year  in  annual  O&M  and  expected  duration  of  5 
years  or  more).  For  sites  with  high  costs  and/or  high  durations,  such  as  a  yet-to-be  constructed 
P&T  system  where  fewer  cost  and  design  parameters  are  fixed,  the  potential  life-cycle  cost 
savings  may  be  even  more  significant.  For  example,  in  the  Blaine  demonstration,  potential  cost 
savings  of  approximately  $10  million  were  identified  relative  to  the  trial-&-error  solutions. 

For  cases  where  the  objective  function  is  not  in  terms  of  life-cycle  cost,  the  cost-benefit 
evaluation  is  less  straightforward  (e.g.,  almost  50%  less  mass  remaining  in  layer  1  for  Umatilla 
using  optimization  algorithms  versus  trial- &- error).  It  is  hard  to  quantify  the  extent  to  which 
benefits  from  reduction  in  mass  remaining,  or  an  increase  in  contaminant  removed,  can  be 
compared  to  the  additional  cost  associated  with  application  of  transport  optimization.  However, 
as  discussed  earlier,  the  premium  of  applying  transport  optimization  algorithms  (up  to  $40,000) 
instead  of  a  trial- &-error  method  is  not  so  high  that  it  would  be  prohibitive  for  most  sites  to 
consider  the  transport  optimization  approach,  and  qualitatively  it  appears  that  use  of  transport 
optimization  should  be  encouraged.  Additional  investment  may  be  required  for  uncertainty 
analysis  of  the  underlying  simulation  model  on  the  optimal  solutions,  but  such  analyses  would  be 
performed  for  the  trial- &- error  method  as  well. 
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6.  Implementation  Issues 


6.1  Cost  Observations 

The  computational  time  for  performing  one  flow  and  transport  simulation  is  critical  during 
transport  optimization  modeling.  Based  on  this  project,  it  is  recommended  that  the  underlying 
flow  and  transport  model  be  modified  to  the  extent  possible  prior  to  optimization  to  shorten  the 
computational  time  without  losing  model  accuracy. 

The  two  transport  optimization  groups  were  from  universities.  It  is  possible  that  the  costs 
incurred  by  universities  might  differ  from  those  associated  with  consultants.  However,  when  the 
optimization  groups  were  selected  a  bid  was  also  received  from  a  consulting  group,  and  the  bid 
costs  were  similar.  Thus,  this  does  not  appear  to  be  significant  issue. 

The  cost  associated  with  the  formulation  process  is  also  a  major  component.  The  formulation 
process  is  inevitable  and  time  consuming  whether  the  trial-&-error  approach  or  the  transport 
optimization  algorithms  are  selected.  It  is  important  that  the  cost  associated  with  the  formulation 
process  be  considered.  In  addition,  the  costs  and  schedule  for  making  any  model  modifications 
in  advance  of  any  optimization  should  also  be  considered. 


6.2  Performance  Observations 

One  issue  in  evaluating  performance  with  respect  to  this  project’s  performance  criteria  was 
“moving  targets”.  That  is,  the  objectives  of  the  installation  or  the  site-specific  transport  models 
changed  by  the  time  the  demonstration  was  completed,  making  the  results  less  relevant  to  the 
installations.  This  was  largely  due  to  the  requirements  for  this  demonstration  project  that  a 
specific  version  of  the  model  be  selected  and  then  used  for  the  duration  of  the  demonstration,  and 
that  the  optimization  modelers  work  independently  from  the  installation  for  an  extended  period 
of  time.  This  issue  points  to  the  benefit  of  closer  collaboration  between  an  optimization  team 
and  an  installation  during  the  optimization  period. 

Minimizing  the  life- cycle  costs  was  always  selected  when  formulating  the  problems  instead  of 
minimizing  annual  O&M  costs  or  minimizing  cleanup  time.  This  is  because  life-cycle  costs  can 
potentially  be  reduced  even  if  annual  O&M  costs  increase  or  cleanup  time  increases.  Therefore, 
minimizing  life-cycle  cost  appears  to  represent  the  more  comprehensive  approach.  For  this 
project,  the  life-cycle  costs  were  calculated  based  on  the  system  duration.  For  two  of  the  cases 
(Umatilla  and  Blaine),  system  duration  was  a  variable  that  was  part  of  the  optimization  problem, 
defined  as  the  time  when  all  concentration  criteria  are  satisfied  in  the  model.  In  reality,  there 
may  be  other  criteria  that  need  to  be  satisfied  (e.g.,  a  certain  period  of  time  where  criteria  remain 
satisfied)  before  a  system  can  be  terminated.  Long-term  monitoring  costs  following  pumping 
shut-off  were  also  not  considered  in  the  life-cycle  cost  evaluations. 
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6.3  Other  Significant  Observations 

The  competitive  nature  of  this  project,  in  which  each  team  worked  in  isolation  from  the 
installation  and  each  other  during  the  entire  modeling  period,  does  not  necessarily  lead  to  the 
best  results  in  terms  of  implementation  feasibility.  For  example,  initial  results  at  the  Umatilla 
site  revealed  that  the  mass  remaining  after  20  years  (according  to  the  model)  would  be  extremely 
low,  so  that  solutions  to  the  third  formulation,  minimizing  mass  remaining  over  20  years,  would 
require  extended  pumping  for  little  or  no  benefit.  Had  interaction  with  the  installation  been 
possible,  this  formulation  would  most  likely  have  been  modified  to  reduce  the  pumping  period. 
These  types  of  issues  reinforce  the  need  for  substantial  interaction  between  the  optimization  team 
and  the  installation  during  application  of  this  technology. 

Another  factor  to  be  considered  during  implementation  is  that  the  optimization  codes  required 
substantial  expertise  to  achieve  successful  implementation  at  the  types  of  complex,  real-world 
sites  considered  in  this  project  (which  are  the  types  of  sites  where  this  technology  is  most 
needed).  This  factor  points  to  a  need  for  technology  transfer  to  train  interested  parties  for 
successful  future  implementation. 


6.4  Lessons  Learned 

The  development  of  mathematical  formulations  of  the  optimization  problems  was  a  difficult  and 
time-consuming  process.  However,  this  formulation  process  results  in  a  concise  and  quantifiable 
statement  of  project  objectives  and  constraints  not  only  necessary  for  transport  optimization 
algorithms  but  useful  for  trial- &- error  method  as  well.  In  that  respect,  the  formulation  process  is 
worthwhile  whether  or  not  mathematical  optimization  algorithms  (or  even  pump  and  treat)  are 
ultimately  applied. 

Some  modifications  to  the  existing  flow  and  transport  model  were  necessary  prior  to 
optimization.  These  modifications  included  changes  to  model  time  discretization  to  correspond 
with  management  periods  in  the  optimization  formulations,  simulating  the  model  under  current 
conditions  into  the  future  to  provide  initial  conditions  for  the  optimization  simulations,  and 
modifying  the  model  solution  package  parameters  to  shorten  computational  time  because  the 
model  runtime  is  the  limiting  factor  for  transport  optimization  algorithms  to  investigate  a  greater 
number  of  potential  solutions. 

Due  to  the  specific  needs  of  this  demonstration  project,  the  optimization  formulations  were  fixed 
at  the  beginning  of  the  simulation  period,  and  simulation  period  length  was  defined.  However, 
normally  the  optimization  modeler  would  interact  with  the  installation  to  develop  revised 
formulations  as  optimization  proceeds,  and  to  adjust  to  new  knowledge.  This  project 
demonstrates  that  such  iterations  should  be  a  useful  component  of  real-world  applications. 
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This  project  also  demonstrated  that  applying  the  transport  optimization  algorithms  was  more  than 
just  “hitting  the  go  button”.  It  required  expertise  to  limit  the  potential  solution  space  to  be 
searched.  If  a  single  optimization  run  were  set  up  to  solve  the  entire  problem  as  formulated,  with 
all  possible  pumping  rates  and  well  locations  in  all  potential  management  periods,  the  number  of 
decision  variables  would  be  much  larger  and  the  computational  times  associated  with  the 
optimization  algorithms  would  be  prohibitive  on  today’s  computers.  Instead,  the  transport 
optimization  teams  employed  sequential  solution  approaches  to  reduce  computational  effort,  in 
which  some  parts  of  the  problem  were  fixed  while  others  were  optimized.  In  some  cases, 
problems  were  solved  one  management  period  at  a  time,  and/or  determining  well  locations  first 
assuming  steady-state  pumping  rates  followed  by  optimizing  well  rates  for  those  pre- determined 
well  locations.  At  some  sites,  surrogate  functions  to  the  simulation  model,  such  as  artificial 
neural  networks,  were  also  used  to  reduce  computational  effort.  The  surrogate  functions,  which 
can  be  evaluated  much  more  quickly  than  the  original  simulation  model,  are  then  used  in  place  of 
the  simulation  model  for  optimization.  All  of  these  approaches  require  expertise  and  professional 
insight  to  be  used  appropriately  without  introducing  significant  error. 


6.5  End-User  Issues 

The  project  teams  for  the  Tooele,  Umatilla,  and  Blaine  systems  were  involved  in  the 
demonstration  project  from  the  screening  phase  through  the  presentation  of  optimization  results. 
These  project  teams  included  the  installation  managers  (for  Umatilla  and  Tooele),  the  USACE 
project  managers  and  technical  staff,  and,  in  the  case  of  the  Blaine  team,  contractors  responsible 
for  the  feasibility  study.  For  the  Umatilla  project,  the  Seattle  District  USACE  Innovative 
Technology  Advocate  was  also  involved.  The  relationships  between  the  ESTCP  team  and  the 
project  teams  were  extremely  constructive.  As  stated  previously,  the  project  teams  provided 
important  input  on  the  formulations  and  were  very  interested  in  the  results.  As  noted  previously, 
in  future  applications  of  the  technology,  closer  contact  with  base  environmental  managers 
throughout  the  optimization  period  would  also  be  essential  to  obtaining  results  that  are  as 
relevant  and  as  up-to-date  as  possible  with  the  installation’s  needs. 

Currently,  both  the  Umatilla  and  Blaine  teams  plan  on  using  the  results  of  the  ESTCP 
demonstration  project  as  a  basis  for  future  operational  changes.  The  Umatilla  project  team  has 
ceased  use  of  one  infiltration  basin  based  on  the  recommendation  of  the  optimization  teams  and 
are  in  the  process  of  seeking  funds  to  update  the  ground  water  flow  and  transport  models  to 
reflect  new  site  characterization  data  before  revisiting  the  optimization  further.  The  Blaine 
project  team  is  considering  the  optimization  recommendations  as  they  proceed  with  preparation 
of  a  Proposed  Plan  and  Record  of  Decision.  They,  too,  are  planning  to  update  their  model  based 
on  recent  minor  site  characterization  efforts.  The  Tooele  project  team  has  been  directed  to 
investigate  the  temporary  (two-year)  termination  of  the  operation  of  the  pump  and  treat  system  to 
evaluate  various  processes  affecting  contaminant  fate.  The  implementation  of  any  of  the 
optimization  recommendations  will  be  postponed  pending  this  evaluation.  Overall,  the 
installations  were  very  open  to  the  recommendations  and  are  implementing  the  recommendations 
to  the  extent  possible  given  other  constraints. 
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This  project  did  not  address  uncertainties  in  the  underlying  groundwater  transport  models.  With 
respect  to  end-users,  it  is  important  to  always  remind  the  end  user  that  the  optimization  results 
are  determined  using  model  simulations,  and  are  subject  to  whatever  uncertainties  may  be 
associated  with  such  predictions.  The  use  of  optimization  algorithms  does  not  increase  or 
decrease  such  uncertainties  relative  to  use  of  trial- &- error. 
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